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ABSTRACT 

It has been a longstanding problem to determine, as far as possible, the characteristic masses of stars 
in terms of fundamental constants; the almost complete invariance of this mass as a function of the 
star-forming environment suggests that this should be possible. Here I provide such a calculation. The 
typical stellar mass is set by the characteristic fragment mass in a star-forming cloud, which depends 
on the cloud's density and temperature structure. Except in the very early universe, the latter is 
determined mainly by the radiation released as matter falls onto seed protostars. The energy yield 
from this process is ultimately set by the properties of deuterium burning in protostellar cores, which 
determines the stars' radii. I show that it is possible to combine these considerations to compute a 
characteristic stellar mass almost entirely in terms of fundamental constants, with an extremely weak 
residual dependence on the interstellar pressure and metallicity. This result not only explains the 
invariance of stellar masses, it resolves a second mystery: why fragmentation of a cold, low-density 
interstellar cloud, a process with no obvious dependence on the properties of nuclear reactions, happens 
to select a stellar mass scale such that stellar cores can ignite hydrogen. Finally, the weak residual 
dependence on the interstellar pressure and metallicity may explain recent observational hints of a 
smaller characteristic mass in the high pressure, high metallicity cores of giant elliptical galaxies. 
Subject headings: ISM: kinematics and dynamics — radiative transfer — stars: formation — stars: 
luminosity function, mass function 



1. INTRODUCTION 

The question of the origin of the stellar mass scale dates 
back to the time of Eddington and Jeans. Modern ob- 
servations reveal that the median stellar mass remains 
unchanged in star forming environments that vary by 
orders of magnitude in density, pressure , metal content, 
and other variables (Bastian ct al. 2010), which suggests 
that this mass must be set mostly by fundamental con- 
stants. A significant advance toward solving the problem 
has come from the realization that the characteristic frag- 
mentation mass of interstellar clouds is d etermined pri- 
marily by their tem perature structures (Sp aans fe Silkl 
120001 : ILarsonl [2005) . In retrospect this is not surpris- 
ing, and can in fact be deduced simply from dimensional 
analysis. An isothermal, turbulent, self-graviting, mag- 
netized gas is fully characterized by three dimensionless 
numbers (e.g. the ratio of gas to magnetic pressure ^ mag , 
the Mach number M., and the virial ratio a v ir)j but 
these ratios admit a rescaling that leaves the dimension- 
less numbers fixed but changes the mass s cale arbitrar- 
ily. ( For a formal proof see the appendix of iMcKee et al.l 
2010). As a result, any simulation or analytic calcu- 
lation of the evolution of an isothermal cloud can al- 
ways be rescaled to produce objects of arbitrary mass, 
and the characteristic mass found in isothermal simula- 
tions depends mos tly on the numerical resolution used 
(|Martel et al.ll20 06). For this reaso n, analytic theories of 
the initial mass function (IMF; e.g . iPadoan fe Nordlundl 
l2002t iHcnnebe lle fe Cha bricr 2008) or simul ations of star 
cluster formation (e.g. lGirichidis et al.ll201 Tl ) that assume 
a purely isothermal equation of state may be able to pre- 
dict the functional form of the stellar mass distribution 
as a function of quantities like A4, /3 m ag, or a V j r , but they 
are always forced to leave the absolute mass scale as a 



free parameter. Thus any explanation of the characteris- 
tic stellar mass and its invariance must somehow depend 
on deviations from isothermality. 

Given this realization, attention has therefore turned 
to the question of what processes determine the gas tem- 
perature structure, and what mass scale they select. 
One idea is that the characteristic stellar mass is set 
by a change from poor t o strong dust-gas coupling at a 
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this explanation faces the problem 
10 4 Mq region of star cluster forma- 
tion in the Mil ky Way is an order of magnitude denser 
than this (e.g. iShirlev et al.l [200l iFaundez et al.ll200l 
iFontani et alll2005D . vet still manages to fragment down 
to masses < Mq . Another proposed mechanism that ap- 
plies to gas near galactic centers or in U LIRGs is heat- 
ing b y either x-rays or cosmic rays (e.g. iKlessen et all 
1 20071: iPapadopoulosI 120101: iHocuk fc Spaansl l2010afl rl 
iMeijerink et al.1 120111 ). While these effects may indeed 
yield deviations from the canonical IMF in extreme con- 
ditions, they do not apply to most star-forming environ- 
ments. Even near the Galactic center, observations so 
far have yet t o find any evidence f or their influence on 
the IMF (e.g. lBrandner et al.ll2008l ). 

Instead, the primary mechanism for determining the 
temperature structure of star-forming gas clouds ap- 
pears to be radiation feedback pro duced by accre- 
tion onto the protostars themselves ()Krumholzl 120061 : 
iWhitehouse fc Ba te 2006). T his is the dominan t energy 
source in a star-forming cloud (lOffner et al.ll2009fl . and at 
mctallicities > 1% of the Solar value ( Omukai et al. 2010; 
Mvers et al1l2011| ) and densities > 10 4 cm" 3 ([Goldsmith! 



2001), which characterize almost all star-forming envi- 
ronments we are able to observe, this energy is well- 
coupled to the gas. Simulations and analytic estimates 
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show that i t changes how gas fragments (iKrumholz et all 
20071 120101 [2011b jKrumholz fc McKed 120081: IBatd 1200a 
Offnc r et al 112009ft . 

The only attempt thus far to understand why radiative 
feedba ck leads to a n invariant peak of the stellar IMF is 
that of Bate (2009), who gives a scaling argument for why 
the density should have little effect on the IMF peak. In 
this paper I e xpand and improve this argument in several 
ways. First, iBatd relies on an empirically-determined 
mass-radius relation to set the protostellar luminosity, 
and it is not clear how this might vary with star-forming 
environment or what physics sets it. I demonstrate that 
the necessary relation can be obtained, at least approx- 
imately, from fundamental constants. Second, the ar- 
gument in IBatd is limited to the case where the gas is 
optically thin. However, the optical depth of a protostel- 
lar core varies radically as a function of wavelength, and 
it is not clear at what wavelength the condition of op- 
tical thinness must be satisfied. Consequently, it is not 
clear to which, if any, star-forming environments Bate's 
argument can be applied. In contrast, the calculation 
I present here requires no assumptions about the opti- 
cal depth of the star-forming region. Third, Bats only 
gives a scaling argument for why the characteristic stel- 
lar mass depends little on the ambient density, but does 
not actually estimate what this mass scale is. I provide 
such an estimate in terms of fundamental constants. The 
argument presented here the refore sign i ficant ly ex pands 
the th eory first advanced by iKrumhoia (|2006l) and IBatd 
(2009) for how radiative feedback can set the stellar mass 
scale. 

2. CALCULATION OF THE CHARACTERISTIC MASS 
2.1. Fragmentation of Interstellar Clouds 

Consider a region of mean pressure P that begins to 
collapse to form a star. Because the dynamical time 
varies with volume density p as id yn ~ l/y/Up, the 
densest portion of a cloud has the shortest dynamical 
time, and rapidly collapses to form a thermal pressure- 
supported seed protostar. The surrounding gas may then 
either accrete onto the existing protostar, or it may col- 
lapse independently to form additional stars. We can 
idealize the fragmentation of this gas as a competition of 
two processes. Near the star, the gas will be heated by 
the star's radiation output, and this will raise its pressure 
and make it resistant to fragmentation. This material is 
therefore likely to accrete. Far from the star, the gas 
is colder, and it is therefore likely to fragment into new 
stars rather than accreting onto the existing one. The 
goal then is to compute the typical mass M of gas that 
is heated to the point where it will accrete rather than 
fragment - in effect to compute the protostar's thermal 
zone of influence. 

In order to make this estimate, I approximate that 
this mass M has a density distribution p w p e {r / R)~ k p , 
where the radius R is to be determined and p e = [(3 — 
k p )/4:Tr]M/R 3 . The choice of k p is arbitrary, and I show 
below that it makes almost no difference. Obviously the 
overall assumption of spherical symmetry is not fully re- 
alistic, but this assumption enables us to perform an an- 
alytic calculation yet still capture the essential physics. 
The virial theorem implies that, for non-magnetized ma- 
terial in virial balance, the mass and radius are related 
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Fig. 1. — Stellar mass M* (blue solid line) and Bonnor-Ebert 
mass times efficiency emMbe (red dashed line) as a function of 
stellar mass M*. The Bonnor-Ebert mass is computed using the 
formalism described by equations J3J - \YI\ . using the fiducial 
parameters k p = 1.5, 5 = 1, j3 = 2, n = 3/2, ejvf = 1/2, 
= 3/4, and computed at a pressure P/kg = 10 7 K cm -3 . 
Note that £m^be 3> for small M*, but that this reverses 
at large M*. The intersection of the two lines gives the esti- 
mated stellar mass set by fragmentation, as indicated by the arrow. 
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where a V i r is the virial ratio. For a marginally bound 
object a v ir ~ 2, and I adopt this value throughout. 
One might worry that the effective pressure might be 
enhanced by the ram pressure of inflow. However, 
observations of the flows around protostars indicates 
that, on the scales of individual low- mass protostellar 
cores, the inflow velocity is at most transonic, indicat- 
ing that the infall ram pr essure cannot be much larger 
than the above estimat e (lAndre et al. 1 12007b iKirk et all 
2007b IRosolowskv et all 120081: iFriesen et all 120091 120101: 



Maruta et all 120101) . This lack or strong shocks within 
cores is consistent with theoretical models of turbu- 
lent f ragmentation (lOffner et al.ll2008b iGong fc Ostrikerl 
2009). One might also worry that, if protostellar cores 
are sufficiently dominated by magnetic pressure, they 
may have a v ir -C 1. In the Milky Way th is does not ap- 
pear to be the case (e.g. lLada et al.l l2008). but we cannot 
directly rule out the possibility that cores are highly mag- 
netic pressure-dominated in other galaxies. One should 
keep this caveat in mind. 

The minimum mass that is capable of undergoing grav- 
itational instability and collapsing to form a second, 
separate protostar is the Bonnor-Ebert mass, A/be = 
1.18Cg/y G 3 p, where c s is the gas sound speed. Nu- 
merical simulations including magnetic fields show that 
they do not significantly alter this characteristic frag- 
ment mass (Pa doan et all [20071 ) . Thus if we consider 
successively larger spheres surrounding the first proto- 
star, the smallest such sphere that contains enough mass 
to be capable of fragmenting and forming another star 
has a mass Mbe- The material interior to that will there- 
fore have to accrete onto the first star (or be ejected 
by its outflow). We may therefore think of the time- 
evolution of the system as follows. The instant after a 
protostar forms, its mass M„ is very small. In contrast, 
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the Bonnor-Ebert mass in the gas around it is 
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and this is much larger than . Consequently, although 
the protostar is small, it has a much larger reservoir of 
gas around it that is too warm to fragment, and will 
instead accrete. As the star gains mass, both and 
Mbe rise, but Mbe rises much more slowly - the star is 
consuming mass faster than its reservoir is growing. Once 
the mass accreted onto the star is equal to the entire mass 
of the heated reservoir, fragmentation becomes likely and 
the star will stop growing. This condition provides our 
estimate for the characteristic stellar mass: 



M* = e M M BE = 1.186am 
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where e M « 1/2 (IMatzner fc McKeel 120001: lAlves et al.l 
120071: Enoch et al.ll2008f ) is the fraction of the mass that 
has collapsed (not the fraction of the entire available mass 
reservoir) that has been incorporated into the star rather 
than being ejected by the protostellar outflow, T e is the 
gas temperature at the edge of the region that will form 
the star, ma is the hydrogen mass, and pa 2 = 2.33 is 
the mean particle mass (in units of toh) for a molecular 
hydrogen-dominated gas of Solar composition. Figure [T] 
illustrates this procedure graphically: at small M*, the 
available reservoir mass cm Mbe M*, but as M* rises, 
eventually the two become equal and fragmentation sets 
in. 

The temperature T e is set by the radiation of the cen- 
tral star. If this star has luminosity L, this i s well- 
approximated by (|Chakrabarti fc McKeel 120051 120081: 
IMvers et al.l[201lh 
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where 7 = 2/? + 4(fc p — 1). Here the dust opacity is taken 
to follow a wavelength-dependence k\ — (5k (A /A)^, 
where 5 is a dimensionless number and we arbitrarily 
dchne k = 0.27 cm 2 g _1 , A = 100 fim, and T = 
hc/\okB = 144 K. With this p arameterization S sa 1 
and 8 ps 2 for Milky Way dust ()Weingartner fc Draind 
l200l . Dust coagulati on can alter 5 by factors of a few 
and 8 by a few tenths (Ossenkopf & Hcnning 1994), but 
I show below that the results are extremely insensitive 
to these variations of this magnitude. The index fcx and 
the dimensionless constant L in turn are given by 
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Note that this calculation assumes that the only signifi- 
cant source of luminosity is the single accreting protostar 
at the center, which is true only if the thermal zones of 
influence of different stars do not ove rlap. Numerical 
simula tions and analytic calculations bv lKrumholz et al.1 
(j 2 1 lh show that this is a good approximation as long 
as the star formation rate in a protocluster where stars 
are forming is < 10% of the mass per free-fall time. All 
observed star-forming regions, b oth within and outside 
the Galaxy, obey this constraint (|Krumholz fc Tanll2007t 
lEvans et al.ll2009l) . The calculation also assumes that the 
heated region is spherically symmetric. Numerical simu- 
lations show that that this is not a bad approximation, 
since the highly diffusive nature of the radiation-matter 
interaction tends to produce fairly round heated regions 
even in the prese nce of asymmetric features suc h as ac- 
cretion disks (e.g. lOffner et^l2^0ll:l5atil200ll) . 

The luminosity of young low mass stars is dominated 
by accretion, so L will be proportional to the accretion 
rate. This is M* ps eMM/td yn , where td ya ~ 1/y/Gp = 
y(3 — k p ) I (3Gp e ) is the dynamical time in the collaps- 
ing region. For this accretion rate, the corresponding 
luminosity is 
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where Af* and i?* arc the stellar mass and radius, ip = 
GMx/R* is the energy yield per unit mass for accreted 
matter, and cl ps 3/4 is the fraction of the accretion 
power t hat goes into light r ather than into driving an 
outflow (jMcKee fc Tanll200i r T l Equations ©, ©, ©, 
and ©, together with ip, fully determine M». Comput- 
ing ip is therefore the next task. 

2.2. The Stellar Mass-Radius Relation from 
Fundamental Physics 

The energy yield from accretion ijj is dictated by a 
number of factors, but the single most important one is 
deuterium burning. For almost all stars this sets during 
accretion, and it forces the stellar core to a nearly fixed 
central temperature. I approximate D-burning stars as 
n = 3/2 polytropes, although different values of n pro- 
duce qualitatively identical results. Assuming the stellar 
core is fully ionized and dominated by ideal gas pressure, 
this implies that tp = T n (kBT c / /L^mn) , where T c is the 
central temperature, fii — 0.61 is the mean mass per par- 
ticle for a fully ionized gas of Solar composition, and T n is 
a dimensionless number that depends on the polytr opic 
index; for n = 3/2, T n = 1.86 (|Chandrasekharill939fh 

I estimate the equilibr ium cen t ral te mperature T c fol- 
lowing the formalism of I Adamsl (|2008[ ). in which a star 
of mass M* and radius i?* is approximated as a poly- 
trope of index n, for which the density distribution fol- 
lows p(£) = p c f n {£), where p c is the central density, 
£ = r/R* is the dimensionless radius, and /(£) is the 

1 Note that cl can also be reduced by episodic accretion that 
delivers some fraction of the final mass in short-duration bursts; 
however, comparisons with the observed protostellar luminosity 
function suggest that the reduction in luminosity during th e non- 
burst phase is modest, only ~ 25% (Offncr & McKcc 2011), and I 
show below that the characteristic mass depends fairly weakly on 
e^. Thus I do not attempt to include this effect. 
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solution to the Lane-Emden equation 
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er = o, 
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with the boundary conditions /(0) = 1, /'(0) = 0. As- 
suming the star is dominated by ideal gas pressure, the 
temperature is then T = T c /(£). From this solution, the 
dimensionless decay rate for the temperature as a func- 
tion of radius is (3 = ^-j" 1 , where £i is defined by the 
condition /(£i) = e _1 ; the dimensionless mass is 



fJ-o 



(10) 



where is defined by the condition /(£*) = 0. 

The next step in the calculation is to estimate the rate 
of nuclear energy gener ation per unit volu me using the 
Laplace approximation (iFowler et al.lll975l) . which gives 



e = Cp 2 & 2 exp(-39), 



(11) 



where C is a constant that depends on the properties 
and abundance of the reactants, 6 = (£ , G/4fcsT) 1 / 3 , 
and Eg = (naZiZ2) 2 1miic 2 is the Gamow energy for 
the reaction. Here Z\ and Z^ are the charges on the 
two reacting nuclei, itir is their reduced mass, and a 
is the fine structure constant. Comp arison with detailed 
calculations of nuclear reaction rates ( Fowler et al.lll975f) 
shows that this is an excellent approximation for the D 
burning reaction with which we are concerned (discussed 
in more detail below) as long as the temperature is <C 
10 9 K. Note that, if we define 9 C = ( J B G /4fc s T c ) 1 / 3 and 
9 = 9 C /(£) _1 / 3 ) the rate of nuclear energy integrated 
over the stellar volume is = CAttR 3 p 2 I(Q c ), where 



/ 2 ™£ 2 9 2 exp(-39) df . 



(12) 



Given these approximations, Adamsl (|2008l ) shows (his 
equation 25) that the conditions of thermal and hydro- 
static balance within the star require that the central 
temperature roughly obey 



/(e c )e- 



45 p Kc CE G h 3 c 2 



Mo 



n + 1 



(13) 

assuming the star is fully ionized so the mean particle 
mass is [limn. Here k c is the opacity at the center of 
the star, and which we approximate as dominated by 
Thompson scattering, so n c = or/wH^l + X)/2, where 
ctt is the Thompson cross section X is the hydrogen mass 
fraction. The results are quite insensitive to changes in k c 
by factors of a few, or even tens. Note that equation (TTBl 
is derived assuming that energy transport is by radiation 
rather than convection, which is true only in part of the 
star during D burning. However, below I compare the 
value of %p derived from this assumption to the results 
of a detailed numerical model that includes convection, 
and show that convection alters if) by at most a factor of 
a few. 

In order to solve equation (TT3)) for the central tem- 
perature 9 C , the final necessary step is to compute the 



10 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

log (M,/M ) 

Fig. 2. — Dimensionless central tempera ture C as a function of 
stellar mass M„ , computed from equation (1 1 3 P . The solid blue line 
shows the fiducial case of (n 
The red dotted line shows n - 
[D/H] = 2 X 1CT 7 . 



= 3/2 polytrope, [D/H] = 2 X 1CT 5 ). 
: 3, and the green dashed line shows 



reaction constant C, defined by 



C 



M ( (AE) 
6 2 exp(-39) \Mi^2W 2 J 



(14) 



where av is the velocity times cross section for the re- 
action, AE is the net energy released per reaction, fix 
and [ii are the mean mass (in units of mn) per reactant 
of species 1 and 2, and the angle brackets indicate aver- 
ages over the Maxwellian velocity distribution of reacting 
particles. Note that, since {av) depends on the tempera- 
ture as O 2 exp(— 39) in the Laplace approximation, this 
quantity is temperature-independent. 

At the ~ 10 6 K temperatures typical of deuterium- 
burning stars, the dominant D burnin g reaction chain by 
a large margin is (IStahler et alJfl980h 



i H ^-2 
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(15) 
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which yields (AE) = 12.6 MeV per D burned, with 
negligible neutrino losses. The first reaction is the 
rate-limiting step. Its Gamow energy is Eq = 
(4/3)7r 2 a 2 TOHC 2 = 0.66 MeV, and the mean masses per 
particle for the two reactant species are /Xi = /ijj = 
mn/X and fi2 = = ^h/[D/H], where [D/H] is the 
abundance ratio of deuterium relative to hydrogen. For 
interstellar gas, an d thus gas at the o nset of D burning, 
[D/H] w 2 x 10" 5 (jStahler et al.lll980l ); the ratio will de- 
cline with time during the main D burning phase, but 
this does not affect 9 C significantly until essentially all 
the D is burned, as I show b elow. For X = 0.71, and the 
value of (av) tabulated by I Fowler et alj (|1975l ). I find 
C = 2.1 X 10 17 ([D/H]/2 x 10" 5 ) cm 5 s" 3 g" 1 . 

Numerically solving equation (|13[) using this value of C 
gives the results shown in Figured] We see that 9 C is al- 
most completely insensitive to changes in the polytropic 
index, and varies by only tens of percent as M» or [D/H] 
vary by orders of magnitude. The central temperature is 
T c k E G /4k B @c- Thus we have 



Tn \ E G 
49 3 / fiim H ' 



(17) 
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Polytrope, n = 3/2 
Polytrope, n = 3 
Numerical 




M/M A 

Fig. 3. — (a) Energy per unit mass t p re leased by accretion onto 
a protostar of mass M, from equation 1171 . for n = 3/2 (solid blue 
line) and n = 3 (red dashed line) polytropes; the solid black line 
shows the results of a detailed st ellar structure calculation (model 
mC5H of Hosokawa ct al. 2011). (b) Ratio of the polytropic 
estimate of ifi to the numerically-determined value, for polytropic 
indices n = 3/2 (solid blue) and n = 3 (dashed red). Note that the 
IHosokawa et all models are initialized to a mass of 0.01 Mq, and 
that the results are highly sensitive to the choice of initial radius 
until the stellar mass reach es several times this value. Thus one 
only should take the Hosokawa ct al. models seriously at masses 
> 0.05 M Q . 



For = M Q , IfindT c = l.Ox 10 6 K, and ip = 2.5xl0 14 
erg g _1 , which agrees to within a factor of ~ 2 with the 
resul t s of detailed stellar str ucture models (jStahler et al.l 
1980; IHosokawa et al.|[201ll) . Figure H shows a more de- 
tailed comparison. As the plot shows, the polytropic 
estimate agrees with the numerical result to better than 
half a dex at all masses > 0.05 Mq, and to better than a 
dex at all masses > 0.01 Mq. I show below that errors of 
this magnitude have very little effect on the final result. 

2.3. The Characteristic Stellar Mass 

Equation (|17p completes the system formed by equa- 
tions ([T]), 0, dl|, and ©, and uniquely specifies M and 
M*. Before proceeding with a numerical solution, how- 
ever, one can gain considerable insight from an approx- 
imate analytic solution. Notice that equations © - (J7J) 
imply that fcx ~ 0.5 and L m 1 as long as 1/2. The 
quantity R represents the ratio of the radius of the dust 
photosphere to the core radius, and this will be < 1 as 
long as the core is opaque enough that stellar photons es- 
cape primarily by diffusing to frequencies where the core 
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Fig. 4. — Characteristic stellar mass M as a function of 
interstellar pressure P, comparing the numerical solution (blue 
solid line) and the analytic approximation (red dashed line; 
equation JT9J). The numerical solution is for k p = 3/2, <5 = 1, 
= 2, n = 3/2, e M = 1/2, and e L = 3/4. The slight upturn 
at very high pressures is associated with the dust photosphere 
moving past the outer edge of the thermal zone of influence around 
the star. Note that the vast majority of star-forming systems in 

the Galaxy lie at P/k B < 10 8 ' 5 K cm 3 , so higher P/kg values 
are realized only in extreme extragalactic environments, if at all. 



optical depth is ~ 1 rather than by diffusing out of the 
core in space while remaining at frequencies where the 
core is optically thick. This is the case for almost all of 
the parameter space relevant to star formation, as can be 
verified readily from a numerical solution I give below. 
For the purposes of analytic approximation, therefore, it 
is reasonable to adopt fcx ~ 0.5 and L ss 1. Similarly, 
O c is nearly independent of the stellar mass. If one takes 
these approximate values as exact, assumes O c is mass- 
independent, and further takes (3 = 2 (as expected for 
most dust models) and k p = 1.5, then it is possible to 
solve the system of equations analytically. After some 
manipulation, the result is 
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where ag = Gm^/hc = 5.91 x 10 is the gravita- 
tional fine structure constant defined for two protons, 
P P = c 7 /hG 2 = 4.63 x 10 114 dyn cm" 2 is the Planck 
pressure, and in the numerical evaluation in equation 
Ip]) I have used e L = 3/4, e M = 1/2, T n = 1.86, and 
9 C = 12.4. Note that the dust abundance S drops out 
of the problem entirely; numerical simulations show th at 
this is an excellent approximation (jMvers et al.ll20lTI ). 

Pressures in star-forming systems cover a 4 — 6 decade 
range, from the relatively diffuse molecular clou ds found 
in nearby dwarf galaxies (jBolatto et al.l |2008|) (surface 
density £ ~ 0.01 g cm -2 , corresponding to a pressure 
P/k B - GT?/k B - 3 x 10 4 K cm" 3 ) to the densest star- 
forming gas clumps seen in the Galaxy (£ ~ 3 g cm~ 2 , 
corresponding to P/ks ~ 3 x 10 9 K cm -3 ). In extra- 
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Fig. 5. — Characteristic stellar mass M as a function of 
interstellar pressure P, for varying parameters. The solid blue line 
is for the same parameters as in Figure [J] The other lines show 
solutions in which one parameter is different: n = 3 (solid purple), 
k p = 1.1 or 2 (dashed and dotted red), S = 0.1 or 0.01 (green 
dashed and dotted), /3 = 1 or 3 (red dashed and dotted), and 

= 3 or 3/16 (equivalent to multiplying ip by 4 or 1/4 relative 
to the fiducial estimate; black dashed and dotted). I use k B = 1.1 
rathe r than k p = 1 because formally the (Chakrabarti & McKce 
2005) approximation becomes singular at k p = 1; however, 
numerical solutions indicate that the results are nearly the 
same as for kp = IT. I do not show the results of varying the 
geometric parameter ejj because this should not vary systemat- 
ically with interstellar environment, and it simply provides an 
overall scaling. Note that, as for Figure [4] the vast majority of 
Galactic star formation occurs at P/k B < 10 8 - 5 K cm 3 . Also 
note that model values that fall below 0.01 Mq should not be 
taken seriously, since this below the estimated mass at which sec - 
ond collapse to stellar density occurs (Masunaga & Inutsuka 2000). 

galacti c stellar systems such as the cores of giant ellip- 
ticals (Ivan Dokkum et al]|2008[ ) and super star clusters 
(jTurner et al.1 12000) we see surface densities that reach 
even higher values of £ ~ 20 g cm -2 . We do not know 
if these systems formed from gas at similarly high sur- 
face densities, but if they did the corresponding pressures 
would be P/ks ~ 10 11 K cm~ 3 . Equation (fT5| predicts 
that even over this very large pressure range, the charac- 
teristic stellar mass should vary by only ~ 1/3 of a dex. 
For comparison, Figure 0] shows an exact numerical solu- 
tion for M* as a function of pressure. We see that, while 
the pressure-dependence is slightly steeper than that pre- 
dicted in the analytic approximation (mainly because fcx 
is not exactly 0.5), the characteristic mass still varies by 
only a decade or so as the pressure varies by more than 
six decades. 

Figure [5] shows numerical solutions for varying stellar 
polytropic indices, density powerlaw indices, metallici- 
ties, dust spectral indices, and radiative energy budgets 
(i.e. values of e_L or Strikingly, most of these factors 
make no significant difference. Varying n, k p , /3, e^, or ip 
within the plausible ranges of variation indicated by the 
upper and lower curves in the Figure produces less than 
a factor of 1.6 change in the characteristic stellar mass 
at all pressures P/ks < 10 10 K cm" 3 . The only factor 
that matters marginally more is metallicity; decreasing 
6 to 0.1, i.e. using a metallicity that is roughly 1/10 the 
Solar value, induces a factor of 2 change in the charac- 

2 The factor of 4 variation is ip or is chosen to encompass the 
error in the value of ip that results from the polytropic approxima- 
tion. This error shown in Figure |3p. 



teristic mass, while using 5 = 0.01 produces a factor of 
3 — 5 shift. This confirms the analytic expectation that 
the properties of the interstellar environment - metallic- 
ity, dust properties, and degree of gas concentration - 
change the characteristic mass very weakly or not at all. 

3. DISCUSSION AND CONCLUSIONS 

The central results of this paper are equation (fl8|) and 
figure |U which describe the characteristic stellar mass in 
terms of the hydrogen mass multiplied by a series of di- 
mensionless factors. Some of these describe the geometry 
of the stellar accretion flow (e_L, cm), the internal struc- 
ture of protostars (T„), and the chemical composition of 
gas (/ih 2 , ^i), and are always ~ 1. Others depend on the 
relative strength of electromagnetic, gravitational, and 
nuclear forces (a, ac); these are fundamental constants. 
The result also depends on 6 C , which describes the en- 
ergy scale in a stellar core in units of the Gamow energy. 
This is set mostly by the properties of the deuterium 
plus hydrogen fusion reaction, which also ultimately de- 
pends on fundamental constants. Finally, the last term 
depends on the interstellar pressure measured in units 
of the Planck pressure; this is the only term that makes 
any reference to interstellar conditions, and there with 
an extraordinarily weak dependence. We can therefore 
understand why the characteristic stellar mass should be 
invariant over such a broad range of conditions: it is 
set almost entirely by fundamental constants, with an 
almost vanishing dependence on interstellar conditions. 

Furthermore, this result naturally explains why the 
stellar mass scale is such that nuclear reactions can be 
ignited in stars. Until deuterium burning begins in stel- 
lar cores, stars contract rapidly as they gain mass, their 
cores heat up, and ip becomes a strongly increasing func- 
tion of mass. During this phase, as stars gain mass their 
thermal zone of influence rapidly expands, since increas- 
ing mass also increases the energy yield from accretion. 
Only once nuclear burning begins and the stellar core 
temperature is stabilized does the energy yield from ac- 
cretion become roughly constant, and the zone of influ- 
ence ceases to expand as rapidly, favoring fragmentation. 
Thus the onset of fragmentation is directly linked to stars 
reaching a mass such that nuclear reactions can begin. 

Finally, I do find a very weak residual dependence of 
the characteristic stellar mass on the interstellar pres- 
sure and metallicity. These effects are small enough that 
they are likely to be masked within a single galaxy, or 
even over a wide range of galaxies of relatively similar 
properties, by random variations in factors like the ac- 
cretion geometry, dust properties, and interstellar pres- 
sures. However, the dependence on pressure and metal- 
licity may produce noticeable variations in samples that 
include galaxies where stars formed under conditions rad- 
ically different than those found today. In particular, I 
find that the characteristic mass decreases weakly but 
noticeably in very high pressure and high metallicity en- 
vironments such as the cores of giant elliptical galaxies. 
There is preliminary evidence for such a bottom-heavy 
IMF based on the presence of unexpectedly strong ab- 
sorption features characteristic of very low mass stars in 
spectra taken from the central portio ns of giant ellipticals 
(|van Dokkum fc Conrovll2010l [20llh . At this point any 
link between this observational result and the theoretical 
one I derive here is necessarily speculative. We have no 
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direct knowledge of the properties of the gas from which 
these stars formed, and it is possible that the pressure 
was less than one would infer from density of the final 
stellar system. Even if the pressures are high, we possess 
of a limited understanding of the physics of star forma- 
tion in such extreme environments. Nonetheless, this 
work points to the need for further investigation of star 
formation at very high pressures, both observationally 
and theoretically. 
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